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Abstract 

In this comprehensible article we develop, following Fantoni and Rosati formalism, a hypernetted 
chain approximation for one dimensional systems of fermions. Our scheme differs from previous 
treatments in the form that the whole set of diagrams is grouped: we do it in terms of non- 
nodal, non-composite and elementary graphs. This choice makes the deduction of equations more 
transparent. Equations for the pair distribution functions of one component systems as well as 
binary mixtures are obtained. We apply they to experimentally realizable quasi-one dimensional 
systems, the so called quantum wires which we model, within Sommerfeld-Pauli spirit, as a ID 
electron gas or as an electron-hole mixture. In order to use our one-dimensional equations we 
consider pair potentials that depend on the wires width. 
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Keywords: Fermi hypernetted chain approximation, quantum wires, Wigner crystallization, coupled quan- 
tum wires. 
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I. INTRODUCTION 



The study of quantum many particle systems amounts to the first half of past century as 
a natural development of quantum mechanics. In general it has evolved along three main 
lines somehow independent one that aims to describe nuclear matter; other addressed to 
the so called quantum liquids ( 3 He and 4 He at low temperatures) and finally a third one 
devoted to the electron gas. This last is a simple model to study several properties of solids 
(Sommerfeld-Pauli model). More specifically, relative to the third line, the properties of the 
ground state of a gas of interacting electrons have been intensively studied under diverse 
approximations[l]. In general these studies were initially oriented towards the 3D electron 
gas. However, the low-dimensional versions of many particles systems in general and of the 
electron gas in particular have been receiving an increasing attention. At the beginning 
the interest was merely academic , but, in last past years, advances in experimental 
research on organic metals 4], carbon nanotubes[5] and semiconductor nanostructures 
have allowed for true realizations of some of these systems giving place to more practical 
motivations to study them. Today, many body systems in low dimension are very commonly 
used to describe real situations in which the movement of the relevant involved particles is 
limited in one, two or even three directions (e.g. quantum wells, quantum wires and quantum 
dots) . 

In this work we focus on one dimensional fermionic systems, particularly the electron 
gas and mixtures electron-hole as models to study, within Sommerfeld-Pauli spirit, some 
of the metallic and semiconductor properties of quantum wires. Our main tool will be 

n 

the irreducible diagrammatic formalism developed by Fantoni and Rosati 8 in order to 

on 

describe 3D many particle systems. There exist in the literature several realizations [9|]-[13( 
of Fantoni-Rosati formalism which differ among them in the way the diagrams are classified 
and grou ped and in the specific form that the complementary energy variational equation 
takes JjJ], [ill, Here we present our own scheme, adapt it to the one dimensional 

case and apply the resulting equations to quantum wires seen as a ID electron gas or as a 
ID mixture of electron and holes in order to study some of their conductor or semiconductor 
properties, respectively. 

In next Section we briefly show the theory on which we base our calculations, It is divided 
into subsections in order to make more clear the presentation. In the first of these subsections 
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we present the model and the background for the diagrammatic expansion of the distribution 
functions within Jastrow approximation for the many particle wave function. Subsection 
B is devoted to to expand the generating function in terms of crude reducible diagrams 
formed by correlation and exchange lines. We also obtain, by functional derivativing the 
generating function, the diagrammatic expansions for the one and two point distribution 
functions. In Subsection C we show how the reducible graphs transforms into irreducible 
ones and classify them in three classes: non-nodal, non-composite and elementary, each one 
being in turn classified according to the type of lines that converge to their root points. 
Two important exact relations among the sums of the diverse class of irreducible diagrams, 
say the relation of van Leeuwen et al. 16] and an Ornstein-Zernike 121] like equation are 
established in subsection D. Additional relationships derived from the composite structure 
of the non-nodal graphs are also presented in this subsection. All these relations, together 
with the variational equation for the energy of Subsection E, define a system of coupled 
equations for the pair correlation function. This system contains the rather difficult to 
calculate elementary diagrams and, according to which of them we include in our calculation, 
we obtain diverse quantum hypernetted chain approximations named QHNC/n, where n is 
the highest order of the elementary diagrams considered. In this work we will restrict to 
the cases n = (that corresponds to neglecting all the elementary diagrams) and n = 4 
that takes into account just the smallest order (elementary diagrams with 4 vertices). A 
scheme to approximate the 4th order elemntary graphs is given in subsection F. Finally, in 
Subsection G we outline the quantum hypernetted equations for binary mixtures of fermions. 

The remainder of the paper (Section III) is used to apply our QHNC equations to de- 
scribe quantum wires as one dimensional systems formed by electrons or by electrons and 
holes. Strictly speaking, quantum wires are not one-dimensional but quasi one-dimensional 
devices, thus in order to use our description we define effective pair potentials that take into 
account the non-zero value of the wire width. Subsection A is reserved to study the micro- 
scopic structure of quantum wires when they are taken as conductors. To this purpose we 
model them by a ID electron gas with effective pair potentials and calculate pair correlation 
functions and structure factors. We observe how these functions change with the density 
and the wire width. In particular we note their variation from a liquid like behavior to the 
one characteristic of a Wigner crystal. In order to check the QHNC results we compare them 
with variational Monte Carlo calculations. The generalization of the QHNC equations to 
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mixtures is applied in Section B to quantum wires seen now as semiconductors, specifically 
as an electron-hole mixture. We take advantage of the fact that the pair correlation between 
an electron and a hole at contact is a measure of the electron-hole recombination rate to 
make contact with photoluminescence experiments. Finally, in Section C, use is made again 
of the QHNC equations for mixtures in order to describe the correlations between the car- 
riers in two parallel coupled quantum wires and show how, in determined cases, an electron 
of one of the wires binds to a hole of the other one so forming a sort of exciton. 



II. THEORY 



A. Pair distribution functions 



We consider a system of N fermions (electrons) moving on a segment of the real axis 
of length L. Actually our interest is in homogeneous, infinite systems at T = in the 
thermodynamic limit (N — ■> oo, L — > oo with N/L — > p, the constant density). If V(xi,Xj) 
is the pair potential between particles, the system Hamiltonian reads 



N 2 N 
i=l i<j 



(1) 



where Pi is the momentum of particle i and m the mass of an electron. Here we restrict ourself 
to the one component case, the generalization to mixtures being sketched in Subsection 
F. The A-particles wave function i/j(xi, xn) verifies Schrodinger equation Hip = Eip. 
For the ideal gas of fermions (for which is V(xi,Xj) = 0) the A-particles wave function 
i>ideai{xi-> ■■■,Xn) is antisymmetric under particles interchange and can be written 



Adeal{xi, -,X N ) = ^2(-l) P P {<Pl(x 1 ).<p2{.X 2 )...(t)N{.X N )} = det [^.(Xj)] 



(2) 



where P is the permutation operator, <f> a (xj) is the one-particle wave function and 
det [4> a . (xj)\ represents the determinant of the Slater matrix 



Mc 



0l(Xi) <j> 1 (x 2 ) ••• <frl( X N) 
<fa(Xl) 02(^2) ••• 4>2{XN) 

4>N(Xl) 4>n(x 2 ) ••• 4>n(xn). 



(3) 
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Frequently, for the one particle functions are used wave planes of the form 

^( Xj ) = 2_ e ^a(j), (4) 



where a(j) denotes the spin contribution and Q the volume that goes to infinity in the 
thermodynamic limit. The allowed momenta fill a segment (Fermi "sphere") of "radius" 
kp = Trp/u, with v — 2 the spin degeneration. 

For systems in which the particles interact (V(xi,Xj) 7^ 0), correlations are induced. A 
convenient manner to handle these correlations is with a trial wave function of the form: 

V>(xi, ...,x N ) = F(xx, ...,x N )4> idea i(xi, ...,x N ), (5) 

The correlation factor F(x\, xn) is assumed to be symmetric under particles permuta- 
tions, so the system statistics is determined just by the anti-symmetric ideal part. An 
appropriate election for F is the Jastrow factorization 

F(x u ...,x N ) = Y[h(xi,Xj) (6) 

i<j 

where the two particle correlation factor f 2 {xi,Xj) goes to zero when the distance X{j = 
smaller than the range of the repulsive part of the pair potential and to one for 
Xij large, denoting the absence of correlations. 

The main objects in our theory are the one and two particles distribution functions: 



p( 1 )( Xl ) = J\r I dX2 ^ -■> dXNll>N(xi, X N )jj N (xi, ...,X N ) ^ 

f dxi,..., dxNi>n(xi, x N )ip N (xi, ...,X N ) 
which normalizes such that J dxip^(xi) = N and 



P {2 \xi x 2 ) = N(N - 1) $ — dxN ^N( x ii ■■■,x N )ip N (x u x N ) 

fdxi,...,dx N iplf{xi,...,x N )'ip N (x 1 ,...,x N ) 



(8) 



The function p^ l '{x\) is the probability density of finding a particle at X\ whereas p^(xi, x 2 ) 
is the probability density of finding two particles at X\ and x 2 , respectively. 
Defining the generating function G {U, /} = log {ip\ if)) as 

N 

G{U, /} = \0g(4) idea l\ Y\ \f(Xi,Xj)\ 2 JJexp[/(Xfc) \lpideal) , (9) 

i<j k=l 
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where U (x) is an auxiliary function which properly chosen simplifies calculations, it can be 
proved that the distribution functions verify 



P 



SG{U,f} 
5U(x k ) 



(10) 




S 2 G{UJ} 



+ pW(x i )pW(x j )-6(x i ,x j )pW(x i ), 



(11) 



P 



SU(xi)SU(xj) 



with 5(xi,Xj) the Dirac delta. 

Directly related to the two particle distribution function is the pair correlation function 
g(xi,X2) defined by: 



For homogeneous systems is p^'(xi) = p and g(xi,x 2 ) = g{xvi)- 
B. Diagrammatic expression for the distribution functions 

In this Subsection we show how the generating function G {U, /} can be expressed in 
terms of graphs whose edges are of two types: ones that come from the Jastrow factors 
(correlation or dynamic bonds) and the others coming from the ideal part of the wave 
function (exchange or statistical bonds). 

Instead of working directly with the Jastrow factors f(xi,Xj) it is convenient to use the 
functions by = b{j(xi,Xj) = f 2 (xi,Xj) — 1 in order to avoid integrals convergence problems. 
This way , the product of correlation factors in Eq. fl5j) can be expressed as products of 
bond functions 



P 



{2) (xi,x 2 ) = p {l \xi)p {1 \x 2 )g{xi,x 2 ). 



(12) 



n i/o*' x i)\ 2 = n^ 1 + h «) = 1 + e 6 « + e e m« + - 



i<j i<j i<j i<j k<l 



1 +x 2 + x 3 + ...+X N 



(13) 



where 
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i<j 



X 3 = ^ ( b ij b jk + bjkhi + b ki b i:j + kjbjkbki) (14) 

i<j<k 



Each term X p can generically be written 

X p = 22 B p (x 1 ,X2,...,x p ) (15) 

i<j<...<p 

with B p (xi, x 2 , ■■■,x p ) being a symmetric function of the coordinates Xi, x 2 , x p . It is 
convenient to write they in second quantization 

X P = ^J dx 1 ...dxpB p (x 1 ,x 2 ,...,Xp)^(x 1 )...^(x p )^(xp)...^(x 1 ) (16) 

where the field operators *ff'(x) and $?(x) create and destroy, respectively, fermions at x. 
Since X p \ipideai) = for p > N, we see that the finite sum is written as an infinite one: 

N N 



(ipideall (1 + X 2 + X 3 + ... + X N ) Y[ exp U (X k ) \lpideal) = (4>ideal\ JJ exp U (x k ) \lp id eal) 

k=l k=l 

1 °° f 

~\- — - ^ ^ / dx\ . . .dx p ]3 p{jic \ , x 2 , • • • , Xp^j x 

N 

x (ip idea i\¥(x 1 )...¥(x p )^(xp)...^(x 1 )Y[expU(x k ) \ip idea i) (17) 

k=l 

The expected values of the field operators can be calculated, using Wick theorem 
in terms of contractions which can be written in the form of matrix elements pij = 

p(U) \xj) of the density matrix operator p(U): 

N 

(4>ideal\ ^{Xj)^(Xi)Y]_expU(x k ) \lpideal) 

= (tffo))- = Pij (18) 



k=l 



(4 'ideal | Yl eX P U ( X k) \ Adeal) 
k=l 

This way we write the argument of the log in Eq. © as: 
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N N 

{^ideal\ Yl \fij\ 2 Yl eX P U ( X k) \^ ideal) = ideal | JJ exp U(x k ) \lp idea l) + 
i<j k=l k=l 

+ (Sum of all labeled Jastrow graphs) (19) 

The parenthesis is the sum of integrals of products of functions b'^s and p'^s. Each of these 
integrals can be associated to a labeled Jastrow graph. 

A labeled Jastrow graph is a set of p points or vertices labeled by coordinates Xi which 
are linked by correlation bonds bij (represented by dashed lines) and/or oriented exchange 
lines pij (here represented by arrows) in such a way that: 

• each vertex is the extreme of at least one correlation bond b^ 

• each pair of vertices is linked at maximum by one bond b^ 

• each vertex has one oriented exchange line pjj arriving to and one leaving from it. 

The contribution of an oriented labeled Jastrow graph is obtained from the following 
rules: 

• each bond b^ that links vertices Xi and Xj contributes with a factor: Xi Xj = 

b(Xi , Xj ) 

• each oriented exchange line leaving vertex Xj and arriving to vertex xi contributes 
with a factor: Xi — > — Xj = (xi \ p(U) \xj) = p^ 

• multiply the product of the above factors by (— l) Up+nL /p\, where n p is the number of 
the exchange lines, til the number of closed loops formed by exchange lines and p is 
the vertices number. 

• integrate over the coordinates that label the vertices. 

Each p-vertices graph belongs to a family of p\ graphs which are obtained, ones from 
the others, by permuting the labeled vertices. All graphs belonging to the same family 
contribute equal to the total sum. In general -because of symmetry- among the p\ graphs 
obtained by permutation of the labeled vertices, there are S identical graphs, so the family 
is formed by just p\/S different graphs. The factor S is called the graphs symmetry number. 
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To avoid considering graphs which give the same contribution, it is customary to just sum 
distinct non-labeled Jastrow graphs multiplied by (— l) np+nL p\/ S . Symbolically:: 

TV N 

{lpideal\ \\ \fij\ 2 \\expU{x k ) \lpideal) = (lpideal\ exp /7 (x fc ) \lpideal) + 
i<j k=l k=l 

+ (Sum of all non-labeled Jastrow graphs) (20) 

Also, in general, a graph T can be decomposed into the product of independent graphs. 
If the decomposition is in va graphs Ta, vb graphs Tb,- ■ ■ , then the complete graph can be 
represented by the product T = -^fy-^fj • • • • The factorials that appear in this last expression 
are symmetry factors that correspond to the interchange of labeled vertices among identical 
connected parts. We obtain the whole sum by adding over all the possible values of z/4, Ub, 
■ ■ ■ . This sum is equal to the exponential of the sum of all distinct connected diagrams. 
Thus, combining this result with Eq.( !20l) we find that the generating function G, defined by 
Eq.(ED reads 

N 

G(U,b) =log(tf |tf> = \og(^ ldeal \ l[expU(x k )\^ deal )+ T k ( 21 ) 

fc=l K=A,B,- 

According to Eqs. fllOl) and (|lip . in order to calculate the one and two particles distribution 
functions from the generating function G(U, b) we must perform the functional derivatives 
with respect to U . Taking into account Eq.(l2Tl). the problem reduces to evaluate derivatives 
of the kind . Since Yk depends on U through a exchange line that contributes with a 
factor (ari| p{U) |x 2 ), we have: 



5T 



K 



5U(x) 



dx\dx2 



5T 



K 



5 {xx\p{U) \x 2 ) 



dx\dx-i 



5{x 1 \p{U) \x 2 ) 5U(x) 

[S(xi, x) (x\ p(U) \x 2 ) - (xi\p(U) \x) (x\ p(U) \x 2 )] 



•5(an|p(/7)|x 2 ) 

The right hand side can be more pictorially written in the form 



(22) 



5T 



K 



6U(x) 



dx 



2x2 




dx\dx 2 x 2 




(23) 



Thus, equations (jTUl) and (TTTj) say that the one and two particle distribution functions 
p l (x) and p 2 (x, y) are equal to the sum of all connected diagrams with one and two external 
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points x and (x, y), respectively. The first diagrams for p 1 (x) and p 2 (x,y) look as: 

,M W = O + O-O + .0-0 + O + • • • (24) 

X X 

o + o-o + <o-o • o • • • • 

X X 

(25) 

C. Reduction to and classification of irreducible graphs 

Among the graphs that appear in Eqs. (1!H1) and fl25l) there are some that have articulation 
points. Removal of an articulation point from a connected graph causes the diagram to 
separate into two or more components, of which at least one contains no root points. A 
graph that is free of articulation points is said to be irreducible. The following generic graph 




has four articulation points labeled a, b, c and d. By successively eliminating the articulation 
points of a given graph, it is always possible to identify a part that contains all the root 
points. We call this part the irreducible part. Any diagram with articulation points can 
be reduced to an irreducible one if, in the irreducible part, we associate certain functions 
to those vertices that were articulation points. These functions "replace" those parts of the 
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graph which were eliminated. For the diagram shown above the irreducible part is: 




y 

The points in this graph which were articulation points are or intersections of correlation 
bonds (like a in the figure) or intersections of exchange lines (like d in the figure). We will 
denote with a black point those which are the intersection of correlation bonds. Given the 
irreducible part of a graph with a black point, we can assume it can be obtained by reducing 
any of the graphs T which can be formed by superimposing to the black point the root point 
of any graph r that has just one root point. The sum of all these graphs T' is, according 
to Eq. [2H p^(x). Thus we replace the sum of all the diagrams T by the irreducible part 
that is common to all of them and whose black point has the density p^\x) associated. 
Analogously, given an irreducible graph which has a field point that were an articulation 
one and which is the intersection of two exchange lines (arrows), we can assume it can be 
obtained by reducing any of the graphs A which can be formed by superimposing to the 
field point the root point of any graph A' that has just one root point but with the condition 
that this point be one of the extremes of a correlation bond. Let A(x) be the sum of all the 



graphs A', it can be seen (c.f. 20] page 85) that 

A(x)=expU int (x)-l (26) 

where Ui nt (x) is the sum of all those graphs that belong to the set whose sum is A(x) but 
which are not products of graphs in the set. It can be demonstrated that the sum of all 
reducible graphs whose irreducible parts differ in the number of interchanges lines of some 
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subdiagram S (as is exemplified in the following figure) 

a a a 




x x 



equals the value of just one irreducible graph (see next figure) where now S corresponds to 
the limiting case of infinite interchange lines that can be though as just a line with associated 
function p^. 




X 



Here pij denotes the elements of matrix p which replace the elements p^ associated to 
each of the previous n interchange lines (n = 3, 4, ...). We have 

p = p + pA — pAp — pApA . . . 
= p + pA — pAp 

= p{U + U mt ) (27) 

From the last equality we see that, by choosing U(x) = —U int (x), the contribution p of the 
interchange lines in the irreducible diagrams equals p(0). Thus, there are no more unknown 
functions and the distribution functions p^(x) and p( 2 \x,y) are formally expressed as the 
sum of all the connected irreducible diagrams without articulation point linked to root points. 
However, although we know how to evaluate any single diagram, since we are treating with 
infinite sums, to obtain p^(x) and p( 2 '(x,y) by simply adding one by one the diagrams 
so calculated is not practical at all. It is convenient to find a way to globally sum all the 
involved diagrams into a closed form. To this end, here we classify the irreducible graphs 
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into three classes. In turn the diagrams belonging to any of these classes, except those that 
we call linear graphs, can be classified according to the kind of lines that arrive or leave 
their root points. It should be mentioned that the linear graphs are connected irreducible 
diagrams that have two root points with just one exchange line leaving from one of them 
and also just one arriving to the other one. All the linear graphs belong to the same group. 

In what follows, we focus on the two point distribution function. Thus all the graphs 
we will be concerned with are Jastrow irreducible graphs with two root points. According 
to Van Leeuwen, Groeneveld y DeBoer 16( we classify these diagrams into nodal, composite 
and elementary graphs with respect to the root points. We say that a diagram is nodal 
(non-nodal) with respect to its root points if it can (can not) be separated into parts with 
the two root points appearing in different components. 

We call a diagram composite (non composite) with respect to a couple of points 1 and 2 
if it is (it is not) composed of two or more parts such that the only points they share are 1 
and 2. 

It is straightforward to see that any nodal diagram is non composite and that every 
composite diagram is non nodal. 

A graph which is both non nodal and non composite with respect to the points 1 and 2 
is said to be elementary with respect to those points. Clearly an elementary diagram must 
have at least four (root plus field) points. It is convenient here to define the order of a graph 
as the total number of points it has. Thus, the smallest order of the elemental graphs is 4. 

According to the kind of lines that arrive to their root point, we will group diagrams into 
three groups that will be indexed (bb), (be) = (eb) and (ee) (b for correlation bonds and e for 
exchange lines) . A root point is called b if only correlation lines converge to it and it is called 
e if there exist an exchange line arriving and leaving the point (to which can additionally 
converge -or not- correlation lines). A graph is of type (i, j) with i,j = b, e if its root points 
are of kind i and j, respectively. 

We denote with Ay the sum of all the diagrams of class A and type (i,j). We will 
use A = N for the non-nodal graphs; A = C for the non-composite and A = E for the 
elementary ones. For example, iV& e represents the sum of all non-nodal diagrams of type 
(be). Thus is convenient to use the matricial notation (A = N, C, E): 
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, Abb(x, y) Ab e (x,y) , 
A(x,y) = | 1 )U ' 1 lU> | (28) 
A eb (x,y) A ee (x,y) 



D. Hypernetted chain equations 

The sum of all connected and irreducible diagrams with root points x y y of kind (ij) 
= 6, e) will be denoted Lij(x,y). The following result of Van Leeuwen, Groeneveld y 



DeBoer 16( relates in matricial form the non-nodal, non-composite and elementary graphs: 

N(x, y) = E(x, y) + b(x, y)I bb + L(x, y) - C(x, y) (29) 

where l bb is the 2x2 matrix which has the element (bb) one and all the other zero. 

All the graphs we call linear belong to just one kind. We denote with l(x,y) the sum of 
all the linear diagrams and with n(x, y), c(x, y) and e(x, y) the sum of all the linear diagrams 
that belong to the non-nodal, non-composite and elementary classes, respectively. They also 
verify the Van Leeuwen, Groeneveld y DeBoer relation: 

n(x,y) = e(x,y) + l(x,y) - c(x,y) (30) 

We show now another relationship among the three classes of irreducible diagrams which, 
together with Eq. (1291) are the basis of a number of liquid state theories: 



L(x, y) = N(x, y) + J dzN{x, z)J{z)L{z, y) (31) 

where 

J(x) 1 J (32) 

V °/ 

In liquid theory Eq. (|3"T1) is usually called the Ornstein-Zernike relation (O-Z) 20|. The 

form of matrix J(x) is due to the fact that two non-nodal diagrams of the type can 

only be linked in series by a field point which has a factor p^'(x) associated. Analogously 
the sum of linear diagrams l(x, y) can be expressed in terms of n(x, y): 

l(x,y) =pM{x,y)+ [l-p {2 \x,y)]n{x,y)[l-l{x,y)} (33) 
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On the other hand, for an homogeneous system so that p^'(x) = p = N/L, Eq. ( 1251) is 
written in terms of Lij = b, e) in the form: 

p {2 \x\ 2 ) - P 2 = P 2 [g{x 12 ) - 1] = p 2 L bb (x 12 ) + pL be {x 12 ) + pL eb (x 12 ) + L ee (xi 2 ) (34) 

If the matrix E of elementary diagrams is known then, replacing into Eq. (l3TT) the non- 
nodal matrix N by the expression given in Eq. (1291) . we have a matricial integral equation 
which would allow to calculate, except for the fact that matrix of non-composite graphs C 
remains unknown, the matrix L and so, according to Eq. (j34p . g(xi 2 ). In order to determine 
the non-composite diagrams we consider the composite structure of the non-nodal diagrams. 
This way we obtain additional relations among non-nodal, non-composite and elementary 
graphs. These relations are not of matricial kind but they depend on the particular type of 
the non-nodal graphs: 



N bb (x n ) = -1 - C bb (x l2 ) + E bb (x 12 ) + b(x 12 ) + f 2 (x 12 ) exp [C bb (x 12 ) - b(x i2 )) (35) 

N be (x 12 ) = E(x 12 ) + L bb (xi 2 )C be (x 12 ) (36) 

A^ ee (xi 2 ) = E ee (xi 2 ) - l(x, y)c(xi 2 )+C be (x 12 )L eb (x 12 ) + C ee (x 12 )L bb (x 12 ) (37) 

n(x 12 ) = e(x 12 ) + l{x 12 )c(x 12 ) (38) 

Using Eqs. flHUj) . flHBl and flHHl) the sums of linear graphs /(xi 2 ), n(xi 2 ) and c(x\ 2 ) are 
eliminated in favor of the sum of the elementary linear diagrams e(xi 2 ) which we are assuming 
as known. 

The relations of Van Leeuwen, Groeneveld y DeBoer ( |29l) . (1301) : the O-Z relations (13T1) . 
( |33l) and the non-nodal structure relations (l35 j) -( l38j) constitute our hypernetted chain equa- 
tions. 

It remains to know the Jastrow factor f 2 {xi 2 ). It can be calculated from the energy 
variation. A convenient way to proceed is to use the hypernetted chain equations in order 
to eliminate f 2 (x\ 2 ) in favor of g(x\ 2 ). We define 

D = L-N 

d = I — n — p, (39) 
so that D and d mean the sum of all the graphs with two root points that have at least one 
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nodal point and use 

g B {x) = 1 + L bb (x) = f(x) exp [E bb (x) + D bb (x)} 
to algebraically eliminate the functions f 2 , N bb y D bb : 



(40) 



N bb (q) = 
D bb {q) = 



g B (x) exp \-E bb (x) - D bb (x)} 

1 2 



+ 



1 - N be (q) 



S(q) l + N ee (q) 



%)" 1 J , 2N ee (q) + 2N be (q) - N&q) 



S(q) 



S(q) + l+N ee (q) 
N ee {q) + N be {q)} fo - N be {q) + N ee {q) 



1+N ee (q) 



(41) 



where S(q) and Nij(q) (ij = bb, be, ee) denote the Fourier transforms of h(x) — g{x) — 1 and 
Nij(x), respectively. 

Therefore our QHNC equations are written in the form: 



D be (q) = -N be (q)+S(q) 
Nl(q) 



1-N be (q) N ee {q) + N be {q) 



D ee (q) 



d(q) 



l + N ee (q) 



+ S(q) 



l + N ee (q) 
N ee (q)+N be (q 



l + N ee (q) 



l+N ee (q) 



—n(q) para \q\ < Icf 

P ara Isl > kp 

N be (x) = g B (x) [E be (x) + D be (x)\ - D be (x) 
N ee (x) = g(x) - D ee (x) 



nix) 



9b{x) 



e(x) + d(x) + 



cr(x 



- d{x) - 



a(x 



v 



g(x) = g B (x) {E ee (x) + D ee (x) + [E be (x) + D be (x)} 2 

-u[e{x) +d(x) +a(x)/u} 2 } 
where v denotes the spin degeneration of particles and a(x) is Slater function. 



(42) 



(43) 
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E. Variational equation for energy 



Eqs. (l4~2l and (l4Uj) is a system of seven equations for the 8 unknowns N be (x), N ee (x), 
n(x), D be (x), D ee (x), d(x), <7b(x) and g(x) so that an additional equation is needed. As such 
we consider the one we obtain from the variation of energy with respect to g(x). Using the 
Jackson-Feenberg formula 2l|, 22| and Eq. (fTTT) together with hypernetted chain equations 
and (I43p . we can write the energy per particle in the form: 



E 

N 



1 



\H\ 



h 2 r 

E F - ^-Po / dx \9b{x) - 1] 
y J dxg(x) 



,<7 



V 



+ 2d(x)V 2 a(x) 



+ 



^ 2 \nf 2 (x) + V(x) 



(44) 



Replacing f 2 (x) in the previous equation by its expression given in Eq. fT4T]) we obtain 



E 

N 



Ep + Eb + £s + £m, 



(45) 



where 



e B = | y dx^fgjx) 



h 2 

—V 2 + 

m 



^6 (g) + 



5(g) - 1 
S(q) 



Es = P j dxEs(x) 

h 2 Po 



dx < [g B (x) - 1] 



V 



+ 2d(x)V 2 a(x) 



g(x) In 



9{x) 
9b(x) 



1 1 _ nV v A [ 2iV ee (g) + 2iV be (g)-iV b 2 e (g) 

n 2. = - a? 2. l 5 ^) " X J j TTaT^ 



2 5(g) 



N ee (q)+N be 


(q) 




2-N 


be (q) + N ee (q) 




1 + iV ee (g) 


2 



(46) 
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The first term has the same form as the energy per particle for a system of bosons. It is 
a functional of the correlation function g(x) and the structure factor S(q). The other two 
terms are specific of fermionic systems and they depend of functions which can be calculated 
from Eqs. (j%2) and (H5j). 

The energy ( )45l) can be taken as functional of just the correlation function g(x) since the 
remainder functions N be (x), N ee (x), d(x), gs(x), etc. can be solved in terms of it by using 
the hypernetted chain equations (j4"2"l) and fj4"3l . Thus g(x) is given by Euler equation: 

HE /N) =0 
8g(x) 

The presence of the functions N be (x), N ee (x), d(x), gB^x) etc. difficult the evaluation of the 
functional derivative. In order to simplify the calculation, we follow a scheme originally pro- 



posed by Lannto and Siemens[12|. To proceed we consider the energy per particle (Eq. 
as a functional of the seven unknowns: E {g(x), D be (x), D ee (x), d(x), N be (x), N ee (x) , n(x)} . 
Then we minimize this functional with respect to the seven independent functions with the 
restriction that hypernetted chain equations f j4*2i) and f j4"3l) hold. This is equivalent to reach 
the extreme of function 



T = E-J2 f ^A,(x) [^(x) - (rfe)^)] - ^ 1 ^ft(g) - M^)] (48) 

i=l J i=l q 

where Xi(x) y fii(q) are Lagrange multipliers, Ni(x) equals N be (x), N ee (x) and n(x) and Di(x) 
equals D be (x), D ee (x) and d(x) for i = 1,2 and 3, respectively. The functions (rhs)i(x) and 
(rhs)^) with % = 1,2 and 3 represent the right hand side of equations f|42|) and (J43]), 
respectively. 

The function T is independently varied with respect to the seven function 
g, D be , D ee , d, N be , N ee y n and with respect to the Lagrange multipliers too. The varia- 
tions with respect to the Lagrange multipliers give the hypernetted equations ( 1421 and ( )43l) 
again. On the other hand, the variations with respect to Di(x) , say with respect to D be (x), 
D ee (x) y d(x), yield: 

dT de s . , ,dNi(x 



+ Hx)w T iTi-Hi{x) = (49) 



dDi(x) dDi(x) lv J dDi(x) 
In Eq. (H9l . the derivatives are calculated using the right hand side of Eqs. (l4"3"j) and 
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Hi(x) is the inverse Fourier transform of /ij(g). The variations with respect to the functions 

Ni(q), say, N be (q), N ee (q) and n(q) give: 



<9T 



dNi(q) dNi(q) 



de M (q) , ~ , ^ dDj(q) ~ 



dNi(q) 



(50) 



The derivatives are calculated by using the right hand side of Equations (1421) and 

Aj(g) is the Fourier transform of Xi(x). Eqs. ( )49l) y ( 1501 are a set of six linear equations for 
the Lagrange multipliers \i(x) y Hi(x). The fact that we have a linear system is what makes 



Lannto and Siemens approach useful 



12|. 



Finally, variation of the functional T with respect to -v/ 'g(x) yields the zero energy scat- 
tering equation: 



h 

V 2 + V(x) + w B (x) + ws(x) + w M (x) 

m 

Here the induced potential wb is the same as for bosons: 



g(x) = 



(51) 



w B {q) 



k 2 q 2 
4m 



E bb (q)+ S(q)-l 



dEM 
dS(q) 



S(q)-l] 2 \2S(q) + 1 



S 2 (q) 



(52) 



For fermions we have the additional induced potentials ws(x) and %(i): 



w M {q) = 2-^— + 2^/i,(g)- F — . (53) 
The derivatives ^f 3 / 3 ? and are evaluated from Eqs. ( 1461) . whereas ^'/^ and ^'/^ 

<99W dS{q) H v ! " <99W aS(q) 

are calculated by using the right hand side of Eqs. ( j4*2l) and ( l4*3l . respectively. 

The procedure to solve the system of equations fj4"2]) . (j4*3l . fjl9|) . fl50|) and fl51]) with Eqs. 
( I5*2"j) and (15*31 is a follows. Firstly we solve Eq. (T5*T1) as for bosons (with = = 

and wb{x) the inverse Fourier transform of Eq. ([5*2"]) ). The correlation function g(x) so 
calculated is used in the hypernetted chain equations ("4*21 and (l43l to obtain by iteration 
the six functions Di and Ni (i = 1, 2, 3). With these functions we solve the linear system of 
equations (I49H50I) for Aj and /ij that allow us to construct the potentials ws(x) and wm{x) 
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(Eq. 15111) . With these potentials we solve again Eq. (1511) and the iteration is continued 
until selfconsistence is attained. In the examples below few iterations (4 or 5) are enough 
to achieve convergence. 

F. Elementary diagrams 

Up to now the matrix E, whose elements are the sum of the different types of elementary 
diagrams, is taken as known. We have already commented about the difficulties involved in 
its calculation which makes this an almost impossible task, so that some kind of approxi- 
mation, that implies to avoid many of the graphs, is necessary. According to the order of 
retained elementary diagrams we have classified the resulting approximation as QHNC/n, 
where n denotes the maximum order. In particular if we ignore all the elementary graphs 
we have the QHNC/0 approximation, which in the literature is often simply called HNC 
In the calculations below, we will take n = 4 that corresponds to the smallest (non-null) 
order. Even for a given order, to sum the whole set of diagrams is very hard, so in our 
calculations we apply the following strategy. For each kind of graph (bb, eb, ee and dd) we 
take as generators the ones shown in Fig. Q] and observe that replacing some of their bonds 
by some of the generalized bonds defined in Fig [2] we recover many of the elementary graphs 
of the same kind and order as the generator one. 




bb eb ee dd 

FIG. 1: Diagrams used as generators of the elementary graphs. 

— =9b(x) - 1 

= g B (x)[E be (x) + D be (x)] 
=>-= = g B (x)[E ee (x) + D ee (x) + (E be (x) + D be (x)) 2 

- u(e(x) + d(x)) 2 } 
-* =g B (x)[e(x) +d(x)] 

FIG. 2: Generalized bonds. 
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Specifically in our calculations in next section we consider the diagrams shown in FigJS] 



E b b,4( X ) = 
E e bA X ) = 



*======# 



+ 



+ 



O O 6 O 6 O 



o o 

W % 6 




+ 



+ 



+ 





+ 



+ 



+ 





+ 



+ 



+ 



6 6 



FIG. 3: Generalized elementary diagrams used in our calculations 



G. QHNC for binary mixtures 

Since several of the applications in Section III involve two component systems, we outline 
here the QHNC equations for binary mixtures. The supra-indices = 1,2) denote 

the species in the mixture. 

The energy equation (Eq. |4~4|) now reads 



E 

N 



+ 2cf (x) W(x 



i J 

£ l - P y j dxgU{z) {"^V 2 [ln^(x) - D*{x)] + V*. 



+ 



(54) 



with /i u the reduced mass. The QHNC equations H2 and H3] yield 
D be (q) = (I + N be (q))(I - N be ) T S(q) [/-(/- N be )(I - N ee )~ 

-N be -i + {i-N ee y l 

D ee (q) = I-N ee -(I + N ee )~ l + S(q) [/-(/- N be )(I + N ee 

- (/ + iVee)-^/ - N be ) T S(q) + (/ + iV ee )- x (/ - iV be ) T %)(/ - N be ){I + iV ei 



—n l (q) para |g| < /c^ 



(»') 2 (g) 



para |g| > kp 
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(55) 



and 



Ki(x) = g%{x) [Dg(x)] - Dg(x) 
NH(x) = g*{x)-D%{x) 



n\x) = g»(x) 



e l {x) + d\x) + 



(Tlx) 



— d\x) - 



(Tlx) 



9 ij (x) = 9%(x) {E»(x) + D%{x) + [El + Z^(x)] 2 - v%(e\x) + d\x) + a*/*/) 2 , 

(56) 

respectively. The symbols N and D denote matrices whose elements correspond to species 
pairs. 

The function whose extreme must be found is 



i=i j,k J 
3 i 

1=1 j,k V 3 « g 



JVf(x) - (r^)f (*) 



(?) - W (9) 



being 



<9T 



dE 



+ Af(x) 



/.f (x) = 



aof(x) aof(x) 1 dDf k (x) 



and 



(57) 



(58) 



£? 9E ~,k () dDt%) 

dNt(q) dNt(q) * dNf^q) 



Af (g) = 



(59) 



the extremal equations. Again, variations of the functional T with respect to y/g*Hx) yields: 



h 



V 2 + V**(x) + w%{x) + w%(x) + w%(x) 



2/i 

where functions w % ^(x),Wg(x) and w l ^(x) are given by 



\f¥Kx) = 



(60) 



H 2 q 2 I 2S u (q) - 3 , S"(q) 2 /m l + S^(q) 2 /m 



+ 



(S U (q)S 22 (q)-S' 2 (q) 2 ) 2 } (61) 
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W 4v^y\ V ij {S n {q)S 22 {q) - S 12 {q) 2 ) 2 ' 



^) = ^ + 2 pf<4fM (as, 

It should be mentioned that a somewhat different version of QHNC for multicomponent 
systems in 3D has been reported by Lantto j^ . 



III. APPLICATIONS TO QUANTUM WIRES 

Given a heterostructure made of several layers of suitable semiconductor materials, by 
means of diverse experimental techniques it has been possible to conveniently modify the 
band structure of some of the layers {(J,!^. The result is that movement of carriers is 



Au 




a b 



FIG. 4: Schematic construction of quantum wires, a) quantum well: the carriers are constrained 
to move in the layer of GaAs. b) Ga ions, which are implanted into the GaAs layer after being 
collimated by a gold mask, modify the layer band structure in such a way that the carriers are 
constrained to move along the quantum wires of width b. 

restricted in one, two or even the three possible directions. In the first case they are confined 
to a slab of a few nanometers width that can be described, within Sommerfeld-Pauli picture, 
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as a quasi-bidimensional system of fermions known in the literature as "quantum well" (see 
Fig. H]a)[25j. If the carriers movement in this quasi-bidimensional structure is even resticted 



in a second direction to a width b t 



ren we have the quasi-onedimensional structures we are 



interested in: the quantum wires 26|], 27| (see Fig. H] b). Typical values of b run from tenths 
to hundreds of nanometers. 

Within this context, our model of quantum wire will be, according to which properties 
one wishes to account for, or a one component system of charged fermions (electrons) or 
a binary mixture whose particles have opposite sign (electrons and holes). These particles 
can move along a quasi-one dimensional region which is practically unlimited in one of the 
directions but transversally very narrow. Here we are thinking in just one wire. However in 
the experimental array there are several wires arranged each one practically parallel to the 
others. In principle it has sense to study an isolated wire ( Fig. [5]) if we assume that the 
interactions among the carriers belonging to different wires are negligible or if we incorporate 
in some way these interactions as a sort of external field acting on the particles of the isolated 
wire. The effect of coupling between two wires is explicitly considered in Subsection C. 



I i i i r i r i i i r i \/ 



FIG. 5: Assuming that the various parallel wires are not coupled, we can isolate one of these 
quasi-unidimensional devices and treat it as a rigorously one dimensional system by using effective 
pair interaction potentials 

Besides, although in principle we can think of an isolated wire as a quasi-one dimensional 
region, actually we wish to describe it from a rigorously one dimensional point of view (Fig. 
[5]). This can be achieved by defining appropriate effective pair potentials. 

In this Section we apply the hypernetted chain equations obtained above to describe 
some properties of quantum wires. In the next Subsection (A) we consider some aspects 
related with their conductor nature, such as the Wigner crystallization. To this we model 
an isolated quantum wire, within Sommerfeld-Pauli picture, as a one dimensional electron 
gas. In Subsection B we study properties that are related with the semiconductor behavior 
instead, particularly the photoluminescence phenomenon. In this case the quantum wire is 
seen as a ID mixture of electrons and holes. Finally in Subsection C we consider a pair of 
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coupled quantum wires in order to study how the correlations of carriers in one quantum 
wire is affected by the presence of a second one. 



A. Quantum wires as conductors 

In this Subsection quantum wires are modeled as a one dimensional electron gas, say 
a rigorously one dimensional system in which the N electrons interact via long range pair 
potentials and freely move along the x axis. For this system we calculate pair correlation 
functions and structure factors in the QHNC approximation and also perform variational 
quantum Monte Carlo (VQMC) in simulations in order to compare our hypernetted chain 
results. 



1. Effective pair potentials 

We consider the effective mass and envelope function approximations and give account 
of the transversal confinement of wires carriers through a parabolic one particle potential of 
the form: 

U(y 1 , Zl ) = #M+J® (64) 

with m* the electron effective mass and b the wire width. 

For this potential the transversal one particle wave functions are those of an harmonic 
oscillator. If the separation between the energy levels is large enough then we can assume 
that the electrons transversally are in the ground state with wave function 

0(!/1 , 2l ) = _L_exp(-*l±£i) (65) 

so that the longitudinal envelope function ip(x) where x = X\ — x 2 must verify the one 
dimensional effective Hamiltonian 

where V e ff is the effective pair potential 

e 2 |0(yi, 2i| 2 |0(y 2 , z 2 \ 2 dy 1 dzidy 2 dz 2 



V eff (x) = -J 



£o\/(x 2 + (y x - y 2 ) 2 + {z x - z 2 ) 2 
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with e being the charge of an electron. 

We see that all the system properties depend on just two parameters: b and r s . The 
parameter r s = (2pa* B )~ 1 (with as = h 2 eo/m*e 2 the Bohr radius for a medium of dielectric 
constant eo) is a measure of the mean distance between two electrons. It is inversely related 
to the density and also measures the system coupling as expressed by the ratio between the 
potential and kinetic energies. 

Eventually we get for V e ff(x) another expression which also captures the essential of the 
interactions between confined electrons and is very frequently used in the literature, the so 
called Schultz potential {28]: 



2. Pair correlation functions 

In the form we have presented the QHNC equations in Section II, the pair correlation 
functions, as defined by EqsJTS] and [HJ are proportional to the probability density of finding 
an electron at a distance x from another one irrespective of their spins. In this subsection 
we compare our QHNC correlation functions with those obtained from VQMC simulations. 
Since is relatively simple to explicitly incorporate the spins in the simulations, we addition- 
ally show the spin dependent pair correlation functions calculated with VQMC so helping to 
better understand the system structural behavior. To define this last correlation functions 
let consider a system of N = iVj + iV| (iVj, number of electrons with spin up and spin 
down, respectively) moving along a segment of length L. Then 



where p s = N s /L is the mean density for electrons with spin s (s =] or !). 

These functions describe the probability density of finding an electron with spin s' sep- 
arated a distance x from another electron with spin s. This way the correlations between 
electrons with parallel and antiparallel spins are distinguished. 

The correlation functions for polarized and non-polarized phases read 




Vx 2 + b 2 ' 



(68) 




(69) 



27 



9n( x ) = w \ 6 ( x l - x \- x )) J 
9ni x ) = w Et^i (%I - ~ a; ))} fully P° larized 




We can combine these functions to define the so called numeric and magnetic correlation 
functions: 



for fully polarized systems. 

Next we show the correlation functions that we have calculated for the effective pair 
potential given by Eq. (I57j) . With respect to the QHNC approximation recall that we have 
mentioned two versions: one in which all the elementary diagrams are ignored. We call it 
QHNC/0. In the other one, denoted QHNC/4, we include the first elementary graphs of 
order 4. In fig [6] these two versions are compared. 

We observe that the inclusion of the elementary diagrams of order 4 becomes the cor- 
relation function more structured. From now on we will take as approximation QHNC the 
version QHNC/4 

As was mentioned, in order to appreciate the goodness of this approximation we compare 
with the corresponding results obtained from VQMC simulations that, in this way, acts as 
a virtual laboratory. 

In figures [7] a and [H] a, we show the results obtained from the QHNC theory and VQMC 
simulations for the correlation function g(x) corresponding to the wire widths 6 = 1 and 4 
(in units of the Bohr radius). We observe that the smaller is the density more structured 
the system is. This effect is more noticeable for small widths. This is reasonable since 
when b increases the system behavior tends to the one of a bidimensional system. The 



g(x) = g nn {x) = - [g n (x) + g n {x)} 
g mm {x) = - \grt(x) -9n( x )\ 



(70) 



for non-polarized systems and 



9nn 



(x) = g mm (x) = g n (x) 



(71) 
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FIG. 6: Comparison of QHNC/0 and QHNC/4 versions for b = 0.1 and a) r s = 1; b) r s = 2. The 
distance x is scaled by the Bohr radius. The pair potential between the electrons is given by Eq. 
(HZ}. 

larger confinement imposed to the electrons in the case b = 1 manifests into a larger Pauli 
repulsion which can clearly be appreciated by comparing the contact values (x = 0). 

The correlation functions g(x) considered in these figures correspond to non polarized 
paramagnetic electron gas. It should be interesting to analyze the correlations between 
pairs of particles with parallel and antiparallel spins. The panel b of figures [7] and [H] show 
g-\l{x), whereas that in the corresponding panel c we see gn(x). These figures were obtained 
just from Monte Carlo simulations because the QHNC approximation, as it is considered 
here, can not allow, in principle, to calculate separately the correlations between electrons 
with parallel spins and those corresponding to electrons with antiparallel spins. In the 
graphics the Pauli exclusion hole for g^(x) is apparent. Besides, from the observation of 
figures b and c, we conclude that the electrons arrange in such a way that spins up alternate 
with spins down. When the density grows (say, when degeneracy grows) the curves g-\i(x) 
show that the system quickly tends to behave as an ideal gas. Obviously, the main difference 
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FIG. 7: Correlation functions of a ID system of electrons interacting via the pair potential given 
by Eq. (f67|) for b = 1 and diverse values of r s : a) averaged spins; 6)parallel spins; c) antiparallel 
spins. 

with g^i(x) is near the origin where the effects of the fermion statistic are more important. 

We have also have calculated the correlations functions for Schultz potential (Eq l68j) In 
figure P we show the results for diverse values of r s and h. It is worth mentioning the 
good agreement of the QHNC and VQMC calculations, better than that achieved with the 
potential given by Eq. ( 1571) 

3. Structure factors and crystallization 

The correlation functions are adequate to describe the structure of the electron gas in 
real space, but essentially are theoretic objects. More directly related to X ray or neutron 
diffraction experiments are the structure factors. These functions are obtained from the 
intensity of the diffracted radiation (previous elimination of the atomic form factor) and 
simply relate to the pair correlation functions via Fourier transforms 
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FIG. 8: Correlation functions of a ID system of electrons interacting via the pair potential given 
by Eq. (|67p for 6 = 4 and diverse values of r s : a) averaged spins; 6)parallel spins; c) antiparallel 
spins. 



1 f dh 

g s Ax) = 1 + -— / —e ikx [SsAq) ~ <U , (72) 

where the subindices denote the spin of the involved particles. 

From the partial structure factors S ss i we define numeric and magnetic structure factors 
which are linear combinations of the spin-up/spin- up and the spin- up/spin-down structure 
factors: 

S nn (k) = S^(k)+S u (k) (73) 
Smm(k) = Stf(k) - S n (k) 

In figures [TU] and [TT] can be seen the structure factors, numeric (panel b) and magnetic 
(panel c), calculated from Monte Carlo simulations, for two values of the wire width and 
diverse densities. In the panel a we show the structure factor (numeric) in the QHNC 
obtained by Fourier transforming g(x). We mention that the lack of the corresponding 
curve for r s = 10 is because convergence in the solution of the QHNC equations was not 
reached in this case. 
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FIG. 9: Correlation functions of a ID system of electrons interacting via the pair potential given 
by Eq. ([68]) for diverse values of r s and b: a) r s = 1, b = 1; b) r s = 1, b = 2; c) r s = 4.5, 6 = 1. 

The most notable feature of curves is the peak located at in the numeric structure 
factor and at 2kp in the magnetic one. This, as we will see below (eqs. 174"! y [75]) . says that, 
for small wire width the charge (or, in this case indistinctly, particle) correlations dominate 
over the spin correlations. At least for the densities we are considering, the effect seems to 
decrease when the width of the wire is increased, that is, when the system is more similar to 
a bidimensional one, and also when the density grows. This behavior is specially apparent 
in fig. [TT] a, where can be clearly observed that for r s = 2 the peak moves at k — 2kp in 
similar way as for electrons moving on a plane( 29|). 

The peak at k — 4kp observed mainly for small widths and densities, can be interpreted as 
indicating the system tendency towards a crystallization of the Wigner type. The wavelength 
corresponding to this value of k is A = p^ 1 = 4, say, the mean distance between particles. 

Wigner crystallization in one dimensional quantum systems with long range interactions 
was firstly studied by Schultz ([28]) by means of bosonization methods( 30], 31|). The main 
conclusion was that long range forces, even if they are weak, cause a state characterized by 
a long range quasi-order more adequate to describe a one dimensional Wigner crystal than 
a liquid. 
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FIG. 10: Structure factors for b = 0.1 and diverse values of r s . a): QHNC; b): VQMC numeric 
;c): VQMC magnetic. Momenta are reduced by the Fermi momentum. 




FIG. 11: Structure factors for b = 4 and diverse values of r s . a): QHNC; b): VQMC numeric ;c): 
VQMC magnetic. Momenta are reduced by the Fermi momentum. 
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The charge-charge correlations as estimated by Schultz is given by: 

(p(x)p(0)) = Ai cos(2k F x f + A 2 cos(4A; F x)e- 4c2v/I ^ + ... (74) 

where Aip are constants that depend on the interaction. The interesting point is the smooth 
variation of the 4kp term, smoother than any power law, showing an incipient charge density 
wave of wave number 4k p- As we have already mentioned, the period of the oscillations, 
4k f, is the mean space between particles, that is, the expected value for a one dimensional 
Wigner crystal. 

On the other hand, the spin-spin correlations are: 

(S(x)S(0)) « B 1 cos(2k F x)- + (75) 

where the lack of the 4kp term should be noted. 

It is worth mentioning that the correlations of charge (eq. 1741) and spin (eq. [TBI) are related 
with the numeric and magnetic structure factors, respectively, via Fourier transforms. 

B. Quantum wires as semiconductors 

In the previous Subsection, quantum wires were taken as just a one dimensional electron 
gas, that is, as conductors seen from the Sommerfeld-Pauli point of view. Here we will 
consider them as semiconductors where the carriers are the electrons of the conduction 
band and the holes left by the electrons that jump from the valence band to the conduction 
one. Recall that in the effective mass approximation the semiconductors carriers effective 
mass depends on the band structure of the considered material and, in general, it will be 
different for electrons and holes. 

We calculate the pair correlation functions using the QHNC equations for binary mix- 
tures of Subsection III Gl In particular, we will consider the electron-hole pair correlations 
at contact and their relationship with the photoluninescence phenomenon experimentally 
observed in quantum wires. Thus the results that we will show can be taken as a true test 
for the QHNC theory against real experiments (not simulations) 

Our system is constituted of electrons and holes with charges e_ = — e + = e, effective 
masses m* and m* h and densities p e = Ph = P- 
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We assume that the particles interact through effective potentials of the form given by 
EqJB7J 



V:j f (x) = eV^exp QJ) erfc (M) ij = e , h (76) 

where, as before, b denotes the wire width. 



1. Contact pair correlation and photoluminescence 

Photoluminescence phenomena are of interest in relation with photodetectors and laser 
devices. Some features of photoluminescence in semiconductors can be explained in terms of 
the radiative recombination of electrons and holes: an electron of the conduction band decays 
into the valence band in indirect form and annihilate with a hole with the correspondent 
emission of radiation. 

The electron-hole recombination rate in a system as the one we are considering here is 
considerably higher than in an hypothetical ideal system in which the carriers do not interact 
among them. This increase in the recombination rate is usually assigned to an enhancement 
factor g e h{0), say the enhancement factor is rightly the electron-hole correlation at contact. 

Since the correlation function at contact g e h{x = 0) is a measure of the probability that 
an electron meets a hole and the photoluminescence intensity is inversely proportional to 
the recombination time r, then the photoluminescence intensity can be studied in terms of 
the contact correlation functions by using the expression (32 ]. 

- = -9eh(0) (77) 

where tq is the radiative recombination rate for the electron-hole pair when interactions are 
turned off. 

In fig. H2] show the QHNC results for the mean life time r of the electron-hole pairs 
as a function of the width h for two quatum wires: GaAs / AlGaAs y Ino^Ga OA7 As/ InP. 
Also the corresponding experimetal points obtained by a group of Sttutgart University are 



shown 
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The effective masses we use are m* = 0.067 m e ; m* h = 0.45 m e for GaAs/ AlGaAs 341 ] 



and m* = 0.041 m e ; m* h = 0.5 m e for Jno.53Gao.47 As/ InP 35j. The one dimensional elec 



tronic density was obtained by fitting an experimental point, the resulting values: n e 
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FIG. 12: Comparison of the curves r vs. b as calculated with QHNC approximation (lines) and 
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those obtained from photoluminescence experiments in quantum wires and reported in ref. 
(symbols). 

1.47xl0 6 [cm -1 ]/ for GaAs / AlGaAs and n e = 1.85 xlO 6 [cm -1 ] for In . 53 Ga 0Ar As/InP, 
are reasonable for these systems. 

We can take the good agreement between theoretical and experimental results as a vali- 
dation of our QHNC version to study experimental realizations of one dimensional fermionic 
systems. It should be mentioned that perturbation calculations within ladder approximation 
previously performed in our group yield similar results [3^], 



C. Coupled quantum wires 

Up to now we have considered just an isolated quantum wire. However we have already 
pointed out that the manufacture process yields a series of nearly parallel wires which are 
separated by nanometric distances. As a first approximation to this scenario we study a 
couple of parallel wires. To this we use again the extension of the QHNC approximation 
we have developed for binary mixtures (Subsection III Gl) . To fix ideas we denote with e the 
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FIG. 13: Pair correlation functions in the QHNC approximation for coupled wires of reduced width 
b = 1.3 at a fixed density (r s = 2) and diverse distance between the wires, a) electron-electron; b) 
electron-hole; c) hole-hole. 

carriers of one of the wires and with h those of the other one. Here we will use Schultz 
potential (eq. l68l) . For fermions of the same wire the effective pair interaction is: 

VllAx) = / } , i = e,h (78) 

e//V ' Vx 2 + b 2 K ' 

whereas between carriers of different wires the interaction will be: 

V%Jx) = - eV i + j = e, h (79) 

e//v ; v^+FT? 

Fig. [13] shows the pair correlation functions inter and intra wires for the same density 
and distinct distances between wires, whereas in figJH] can be observed the same functions 
but for different densities when the distance remains constant. In these figures the distance 
are in units of the inverse of the Fermi momentum. The ratio of the effective masses to the 
electron mass are = 7 and wire diameter in units of the Bohr radius an = -^f% is b/an 
= 1.3. We are taking a dielectric constant e = 13. In relation with these two figures it must 
be understood that for the functions g e h{x) (panel b) the origin corresponds to a distance d 
(minimum distance between wires). 
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FIG. 14: Pair correlation functions in the QHNC approximation for coupled wires of reduced width 
b = 1.3 at a fixed distance between the wires {d = 1) and diverse densities, a) electron-electron; 
b) electron-hole; c) hole-hole. 

From fig. [13] a y b we conclude that, at least for the values we are considering here, the 
variation in the distance between wires hardly affects the correlations inside a given wire. 
As far as g e h{x) (panel b) its variation with d is more noticeable. This is reasonable because 
of the dependence of the effective pair potential with the distance between wires (see eq. 

Also from fig. [THwe see that, for a fixed distance between wires, the correlations inside a 
given wire (panels a and c) depend on the density as it would be expected; the greater the 
density the smaller the structure as for a degenerate gas. On the other hand the curves for 
g e h(x) show as remarkable feature a noticeable increasing at contact (x = d) when density 
decreases. This fact can be interpreted as the indication of an exciton formation. 

Similar QHNC calculations but for coup led electron and hole quantum wells have been 



reported by Alatalo et. al. [38|], 39|, 40j. More specifically, Yurtsever and Tanatar 41] 
studied coupled quantum wire systems using a perturbation approach based in the ladder 
approximation of Yasuhara [42] . 
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IV. CONCLUSIONS 



In this work we have tackled the problem of describing one dimensional systems of many 
fermions. To this end we have developed a realization of Fantoni-Rosati formalism, which 
yields a version of the hypernetted chain approximation that shows, with respect to other 
versions available in the literature, some remarkable aspects, particularly the way in which 
the distinct classes of graphs are ordered and summed and also the form considered for 
the energy variational equation. We apply our QHNC equations to study quantum wires 
modeled, according to Sommerfeld-Pauli picture, as a one dimensional electron gas or as 
a mixture of electrons and holes in ID. The confinement effects are taken into account 
through one dimensional pair potentials that include as a parameter the wire width. Our 
QHNC results were compared with variational Monte Carlo simulations and, in the case of 
electron-hole mixtures, also with photoluminescence experiments in quantum wires. 
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